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General Properties of Response Functions of Nonequilibrium Steady States 
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We derive general properties, which hold for both quantum and classical systems, of re- 
sponse functions of nonequilibrium steady states. We clarify differences from those of equi- 
librium states. In particular, sum rules and asymptotic behaviors are derived, and their 
implications are discussed. Since almost no assumptions are made, our results are applica- 
ble to diverse physical systems. We also demonstrate our results by a molecular dynamics 
simulation of a many-body interacting system. 
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Stable states of macroscopic systems can be well characterized by their responses to ex- 
ternal probe fields. When the probe fields are weak, the responses are linear functions of the 
probe fields. General properties of the linear response functions are well-known for equilibrium 
states.^' ^ In contrast, those for nonequilibrium steady states (NESSs) are not well understood 
yet, although many attempts have been made.^"^^ 

For example, perturbation expansion of the density operator pp oi NESSs in terms of a 
driving field F has often been employed.^ This gives linear and higher-order (n = 2, 3, • • • ) 
response functions (denoted by <I>cq and <I>eq^ , respectively) of equilibrium states. General prop- 
erties, such as symmetries, of $eq^ were thus derived.^ Similar results were also obtained from 
the fluctuation theorems.^ However, it is generally hard to obtain linear response functions 
$i? of NESSs from <I>cq and $eq^ because such expansion converges only slowly (or does not 
converge) for F that is large enough to drive NESSs of interest. Another approach is to utilize 
some general expression of ^p. Such an expression was derived, e.g., in ref. 7. However, it 
contains the expectation value of a function which is unknown except for simple cases. To 
derive physical results for from such a formal expression, simplifying assumptions were 
made,'' at the expense of generality. Furthermore, one may expect that ^p could be expressed 
by small fluctuation in NESSs. However, refs. 8 and 9 showed that ^p of finite macroscopic 
systems is not a universal function of the fluctuation and temperature, i.e., ^p depends also 
on another system-dependent parameter(s). Because of these difficulties, general properties of 
^p were not clarified. 
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In this paper, wc derive general properties of ^p, which hold for diverse physical systems. 
We clarify which properties are common or different between $f and $cq- We also illustrate 
some of the properties by a molecular dynamics (MD) simulation of a many-body system. 

Response function of NESS - Suppose that a strong static field F is applied to the target 
system (the macroscopic system of interest), and a NESS is realized for a sufficiently long time, 
i.e., for [tinjiout]) where tout — ^in is macroscopically long. In such a NESS every macroscopic 
variable A takes a constant value {A) p in the sense that its expectation value at time t behaves 
as 

{AYp = {A)p + o{{A),^). (1) 

Here, (^)tp denotes a typical value of A, and o((A)tp) represents a (time-dependent) term 
which is negligibly small in the sense that o((^)tp) /(^)tp — >■ as F — >■ oo, where V denotes 
the volume of the target system. When A is the energy [/, for example, (i7)tp = OiV) and 
{UYp = {U)f + o{V). 

Suppose that a weak and time-dependent prohe field fit) is applied, in addition to F, to 
the target system for t > to, where t-m < to < t < tout- We are interested in the response of 
the NESS to f{t). Specifically, we focus on the response, 

AAit) ^ {A)'^^f - {A)p, (2) 

of a macroscopic variable A of the target system. To the linear order in /, AA{t) can be 
expressed as 

AA{t)= f^F{t-t')f{t')dt', (3) 

Jto 

which we call the linear response relation. This and the causality relation 

$F(r) = for r < 0, (4) 

define the response function $F(r) of the NESS, as in the case of ^eq{T).^'^ 

Microscopic expression of^p- Equations (3) and (4) do not refer to microscopic physics 
at all - they are phenomenological equations which are closed in a macroscopic level. We now 
relate them to microscopic physics by deriving a microscopic expression of ^p{t). 

Since we are interested in general properties of NESSs, we do not employ perturbation 
expansion with respect to F,^'^ which, for large of interest, converges only slowly or does 
not converge except in limited physical situations. To treat F non-perturbatively, we consider 
a large system which includes the target system, a driving source that generates F, and a 
heat reservoir (s). We call this large system the total system, and denote its Hamiltonian by 
H^°^. When the target system is an electrical conductor, for example, the driving source may 
be a battery, the heat reservoir may be the air, and the total system is the one that includes 
them all, as shown in Fig. 1. On the other hand, we do not include the source of the probe 
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Fig. 1. An example of a large system which we call the total system. It includes an electrical con- 
ductor, a battery, a heat reservoir, and so on. 



field /, such as a microwave generator, in the total system. We assume that / gives rise to 
the interaction term —Bf{t), where B is a macroscopic variable of the target system. Hence, 
the total system is an isolated system except that it is subject to an external weak field /. 
Therefore, the density operator of the total system evolves as 

^h^/^XJit) = [h'"' - Bf{t), p'^lfit)] . (5) 

We denote p^pXj{t) with / = by p^p^{t). When / = 0, the reduced density operator of 
the target system is 

pp^Tr'[p'^\t)], (6) 

where Tr' denotes the trace operation over the degrees of freedom other than those of the 
target system. For a time interval [tin, ^out], a NESS is realized in the target system (while the 
driving source such as a battery is not in a steady state), and, according to eq. (1), we can 
regard pp as being independent of t as far as macroscopic variables are concerned. Unlike the 
equilibrium case, however, an explicit form of pp is unknown. 

When / 7^ 0, Pp^{t) is changed into p^p^ji^t). We assume that the NESS is stable against 
small perturbations. That is, after small perturbations are removed the target system returns 
to the same NESS as that before they were applied. Except for NESSs of some soft matters 
and glass near a melting point, most NESSs, including those of nonlinear optical materials^' ^ 
and electrical conductors, satisfy this assumption. 

For such NESSs, we can evaluate <5f by evaluating the solution of eq. (5) using a first-order 
perturbation expansion with respect to /. For an observable of interest. A, of the target system, 
its response AA{t) = Tr[p*?^^(t)i] - Tr[y5*/(to)i] = Tr[p*?^_^(t)i] - Tt:[pfA] is evaluated to 
the linear order in / as 

AA{t) = I' ^Tr (^p'^\t') [b, A{t - t')] ) f{t')dt', (7) 
where the symbol ' " ' denotes an operator in the interaction picture, i.e., A{t) = 
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efi 1 . Prom consistency with the macroscopic physics, eq. (3), t m Pp [t ) m 

eq. (7) must be irrelevant. Hence, we can take t' to be an arbitrary time (such as to) in 
[iin,iout]) and simply write p^p^{t') as p^*. We thus obtain a general formula; 

$^S(r) = ^Tr f p*/ [5, i(r)l ) for r > 0. (8) 
in \ I J/ 

Here, we denote by $p^(r) to designate variables A and B. For classical systems, 

[B, A{T)]/ih (and similar expressions in the following equations) should be replaced with the 
corresponding Poisson bracket. 

The right-hand side (rhs) of eq. (8) represents some correlation in the NESS. If it could 
reduce to the symmetrized time correlation,^ it would be equivalent to fluctuation (in the 
classical regime, ksT » hw^^). However, as will be discussed shortly, this is not the case 
when F ^ 0. Hence, we do not call eq. (8) a fluctuation-dissipation relation (FDR). We 
call it the response-correlation relation (RCR). Among similar formulas, eq. (8) has the most 
convenient form to derive useful properties which will be presented below. 

Fluctuation-dissipation and reciprocal relations - Before deriving useful properties, we 
use the RCR to discuss why the FDR and the reciprocal relations (including those for finite 
frequencies^) are violated in NESSs.^'^ 

When F = 0, and ^p^ of eq. (8) reduce to the equilibrium state peq and its response 
function , respectively, and the RCR reduces to the equilibrium one^'^ (which is not 
customarily called the RCR, though). When the canonical ensemble, e.g., is employed, both 
Peq oc e~''^^*°* and e^^*°^^ (which defines A{t)) are exponential functions of H^°^. Using this 
fact, we can recast the equilibrium RCR as 

<''(r) = ^(^(0);i(T))eq, (9) 

where ( • ; • )eq denotes the canonical correlation.^ This result is known as the Kubo formula, 
from which one can derive the reciprocal relations.^' ^ In the classical regime {ksT ^ Hlo), the 
canonical correlation reduces to the symmetrized time correlation,^ and hence to fluctuation,^^ 
and one obtains the FDR.^'^ 

When F ^ 0, in contrast, p^* (although its explicit form is unknown) cannot be an 
exponential function of H^°^ only. As a result, the RCR cannot be rewritten into a form 
similar to eq. (9). That is, the RCR holds both for equilibrium states and for NESSs, whereas 
it is equivalent to the Kubo formula only for the former. As a result, the FDR and the 
reciprocal relations are violated in NESSs. The difference between the rhs of eq. (8) and the 
symmetrized time correlation divided by ksT is the violating term. 

Properties derived from the phenomenological equations - Equations (3) and (4) take the 
same forms as those for <I>cq- Therefore, among many properties of <I>cq, those derivable only 
from eqs. (3) and (4) (without using the Kubo formula) hold also for ^p. For completeness, 
we mention such properties, although some of them may be rather obvious. 
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For stable NESSs, we expect that |$i?(T)| is integrable over (—00, +00).^^ Hence, the 
Fourier transform 



/oo roo 
-00 Jo 



(10) 



should be a continuous function of u. As in the case of the Fourier transform E;eq(a;) of $eq(i" 
we can easily show that "Bpioj) satisfies the dispersion relations, 
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and the moment sum rules. -"^ We also see that ReEpiuj) is even, whereas ImSi^(a;) is odd. 

Properties derived from the RCR - We now present the most important results of this 
paper. For F = 0, many properties were previously derived for $eq from the Kubo formula. ^'^ 
As discussed above, some of them (such as the FDR) are violated for $f when F j^O. However, 
the other properties of $eq can actually be derived from the RCR without using the Kubo 
formula, although they were often derived from the Kubo formula (or similar expressions) in 
the literature. Such properties hold also for ^i;' if some quantities are replaced with those of 
a NESS (see below) because the RCR holds even when F ^ 0. We now present them. 

Note that although their forms are similar to those for $eq/'^ their values are often 
different from those of $eq) as will be illustrated later. 

Integration of (cj) yields 



J — ( 



E^^'ico)^ = $^^(+0) = ^Tr (p*/ 



dw 



B 



,i(0)]). 



Since A(Q) (= A) and B are operators of the target system. 



Tr p'/ S,A(0) 



TlipF 



B,A 



B,A 



where {■)f ^ Tr(P-F " ) denotes the expectation value in the NESS. Noting also that 
Im 'Bp^ [u) is an odd function, we obtain the following sum rule for Re S 



F ' 



J — ( 



ReS^^(a;) 



duj 



B,A 



(13) 



1 

TT \ih 

Since the rhs is the expectation value of a known operator, it can easily be measured experi- 
mentally}'^ 

Moreover, by integrating eq. (10) by parts, multiplying the result with ui, and integrating 
the resultant equation, we obtain the following sum rule for ImH^^; 



L^ImS^^(a;)-( ^ 



1 



ih - 



B,A 
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TT 
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(14) 
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In the second line, we have replaced /5^"* with pp because ^(0) (= [^4, H^°^]/ih) is localized in 
the target system. (Recall that all physical interactions in H^°^ should be local interactions.) 
The expectation values in this sum rule can also be measured experimentally. 

From these sum rules we can see the asymptotic behavior of 'r.p^{uj). As \uj\ is increased, 
^.^^{uj) should decay quickly enough such that the integrals of eqs. (13) and (14) converge. 
In particular, we find 

lim ujImEi^(u) = {^\b,A]) . (15) 

oj^oo ^ ^ ' \in I i / p 

As discussed above, the reciprocal relation for ^^^[t) does not hold for ^p^(t), i.e., for 
for eachuj. However, eq. (13) yields 

/oo poo 
Re^i^ {u)duj = - ReS|^(w)dw, (16) 
-oo J —oo 

i.e., a reciprocal relation holds for the integrated values. 

Implications - The left-hand side of eq. (13) equals to $^^(+0) (representing instanta- 
neous response), which however is hardly measurable in real physical systems. In contrast, 
He.'Ep^ioj) is measurable in a certain finite range of u. For higher which is out of such a 
range, Re!E!p^(a;) decays quickly, as mentioned above. Hence, one does not necessarily have to 
measure it for higher uj. Therefore, eq. (13) should be considered as a prediction on ReE!^^(a;) 
in a certain finite range of uj. For similar reasons, eqs. (14)-(16) become important when one 
wants to get information on 'E^^iu). 

These equations, like the corresponding ones for S^^(a;), are very useful for measuring 
or theoretically calculating H^^(a;). For example, one can check experimental or theoretical 
results against them. Using them and eqs. (11) and (12), one can also estimate S^^iuj) in some 
range of uj from existing data of S^^iuj) in another range. Moreover, as will be illustrated 
for a Langevin model later, our results can show that some equality is identical to another 
equality, which were previously treated as independent equalities. 

Furthermore, we can see the following. According to eq. (13), the sum value (integral) of 
Re Si^(ci;) equals to the expectation value (C) p of the Hermitian operator C =[B, A]/ih. Since 
this is an equal-time commutator, C depends neither on the Hamiltonian nor on the state. 
There is no difference in C between free particles and interacting ones or between equilibrium 
states and NESSs. Only through pp the sum value can be affected by these factors. 

When A and B are linear functions of canonical variables, in particular, (7 oc 1 (identity 
operator) and hence the sum takes the same value for every state. More generally, we can say 
the same when Tr[p^p^{t)C] is conserved during evolution from an equilibrium sate to NESSs 
of interest. 

For example, suppose that the target system is an electrical conductor of length L. A static 
electric field F is applied in the x direction (along the conductor). Let and pi be the x 
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components of the position and momentum, respectively, of the jth electron in the conductor. 
Then, B = eqi, and the electric current averaged over the x direction may be given by 

e 



j 

where m is electron's mass. If one is interested in /, putting A = I yields C = (e^Ne/mL)!, 
where Np is the number of electrons in the conductor. We thus find 



J — ( 



ReS-M^ = ^, (18) 



which is independent of -F. Generally, Re E^^^uj) at low u depends strongly on F for large |-F|. 
At high oj, on the other hand, it is expected that Re E^^{uj) would be insensitive to F because 
each particle would not collide with other particles in a short time period ~ 1/uj. (These 
facts will be illustrated later.) Prom these viewpoints, eq. (18) may be counterintuitive, and 
therefore is useful. 

Regarding ImS^^, we can apply eqs. (14) and (15). For example, the latter yields 

g2 7y 

lim ojImEiPiu) = — (19) 

which is also independent of F. 

For more general cases where Tr[/}^*(t)C] is not conserved during evolution from an equi- 
librium sate to NESSs of interest, the sum value generally depends on F. For example, if 
one is interested in F (to investigate, e.g., current fluctuation) in the above example, putting 
A = P yields C = {2e^Ne/mL)i. We thus find 



" ReS^'^(a;)— = lim a;ImS^'^(a;) = ^^{I)f, (20) 

^ TT oj^oo mL 

which depends strongly on F. This fact demonstrates that although the forms of eqs. (13) 
and (15) are similar to the corresponding ones for ^eq;^'^ their values can be very different. 

Non-Hamiltonian systems - We have assumed that the total system, such as Fig. 1, 
is a Hamiltonian system. In studies of NESSs, non-Hamiltonian models, such as stochastic 
models, are often employed. The general properties of Ep^ must hold also in such models if 
the models are physically reasonable ones, because every existing physical system is believed 
to be a Hamiltonian system if a sufficiently large system (such as Fig. 1) is considered. 

For example, a nonlinear Langevin model 

p = -{j/m)p - U'{x) +F + fit) + m, (21) 

where p = mx and U{x) is a potential, may be derived from a Hamiltonian model by making 
projection and by approximating 7 and ^(t) as a constant and white noise, respectively. If 
these approximations are physically reasonable, eq. (13) must hold in this Langevin model 
because it holds in the original Hamiltonian model. We can prove that this is the case for 
any value of F. Hence, the Langevin model is physically reasonable in view of eq. (13). On 
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the other hand, it is wcU-known that the Langevin model gets worse as \F\ is increased. This 
ihustrates that eq. (13), and the other general properties derived above, are not sufficient hut 
necessary conditions for good nonequilibrium models. In this respect, they are similar to the 
charge conservation, which is also a universal and necessary condition for good models. 

As an illustration of significance of eq. (13) on the Langevin model, we can show using 
eq. (13) that equality (5) of ref. 10 on dissipation {J)f is identical to a simple relation, 
{J)p = {-f/m){{p^)F/m - kBT), of refs. 17 [eq. (82)] and 18 [eq. (4.13)]. 

Numerical example - Finally, we illustrate the validity of eqs. (18) and (19) by an MD 
simulation of a model of a classical two-dimensional electrical conductor.^' ^^"^^ The model 
includes A'^e electrons (e), Np phonons (p) and Ni impurities (i), where Ne,Np,Ni ^ 1. The 
e-e, e-p, e-i, p-p and p-i interactions are all present, whereas the static electric field F acts 
only on electrons. The energy of this many-body system is dissipated through thermal walls 
(which simulate a heat reservoir) for phonons, and a NESS is realized for each value of F. 
This model is a mechanical model supplemented by the thermal walls for phonons. For NESSs 
at large |-F|, we previously found the following: (i) {I)f is nonlinear in F,^^ (ii) the long-time 
tail is strongly modified, ^^'"'^'^ and (iii) the FDR is significantly violated.^ 

To compute H^^(a;), we take the probe field f{t) to be monochromatic; f{t) oc sin(a;t). 
We apply it in the x direction in addition to F, and perform an MD simulation, in which we 
here take e = m = l,A''e = Np = 1500, iVj = 500, and L = 750. By calculating the classical 
counterpart of I{t) of eq. (17) for a sufficiently long time (^ l/i^)) we obtain S^^(w). This 
procedure is repeated for various values of oj and F. 

Figure 2 shows a;-dependencc of Re Sj,^(c<^')/7r for three difi'crcnt values of F. At low 
frequencies ReH^^(a;) depends strongly on F, implying that the response to F is nonlinear. 
At high frequencies, the F-dependence looks quite weak. However, since the horizontal axis 
is in the logarithmic scale and f 'B{u})du! = f H(a;)(i(lna;), small differences in ReH^^(a;) at 
high frequencies contribute significantly to the tj-integral of eq. (18). As a result, the integral 
over all uj (i.e., the left-hand side of eq. (18)) becomes independent of F (within possible 
numerical errors), as shown in the inset of Fig. 2. Moreover, the value of the integral agrees 
with that predicted by eq. (18). 

Furthermore, Fig. 3 shows w-dependence of a;ImS^^(c<;) for three different values of F. 
We observe that as a; is increased oolmE^^ (uj) approaches the same asymptotic value, which 
agrees with that predicted by eq. (19), for all values of F. 

We have thus confirmed eqs. (18) and (19), which have been derived from the general 
results, eqs. (13) and (15), respectively. Conversely, as discussed above, the agreement of our 
numerical results with the general results indicates the following: (i) this model is physically 
reasonable, and (ii) our MD simulation well describes NESSs and their responses. That is, 
validity of the numerical results have been checked against the general results even for large 
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Fig. 2. ReS^^(a;)/7r for F = (circles), 0.06 (squares) and 0.1 (triangles). The inset shows the 
w-integral of ReH^^(a;)/7r plotted against F, where the dashed line represents the rhs of eq. (18). 



In contrast, a typical conventional method is to check results against the FDR (see ref. 11 
and references cited therein), which holds only for small \F\. 

Summary - We have derived general properties of linear response functions $ f (and their 
Fourier transform 'Bp) of NESSs, which are driven by a strong field F. For completeness, we 
have presented all the basic properties (including rather obvious ones, such as the dispersion 
relations) which correspond to those of response functions $eq (^eq) of equilibrium states.^' ^ 
Specifically, we have derived the response-correlation relation, eq. (8), which however cannot 
be recast into the form of the Kubo formula when F 7^ 0. As a result, the FDR and reciprocal 
relations are violated in NESSs, although the latter holds for the integrated values, eq. (16). 
In contrast, the dispersion relations, eqs. (11) and (12), and the moment sum rules hold 
even when F because they come from the phenomenological equations, eqs. (3) and (4). 
Furthermore, the sum rules and asymptotic behaviors, eqs. (13)-(15), hold even when F ^ ii 
the expectation values in an equilibrium state, Tr (peq • ), are replaced with those in a NESS, 
Tr {pp ■ ). We have illustrated some of these results by an MD simulation of an electrical 
conductor. 

These results are quite general, which apply to diverse physical systems because no as- 
sumption has been made except that the NESSs are stable. Further generalization to the case 
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Fig. 3. CO ImS^^(cij) for F = (circles), 0.06 (squares) and 0.1 (triangles). The dashed line represents 
the rhs of eq. (19). 



of time-dependent and/or spatially- varying F (and/or /) is straightforward. 
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